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Some general properties of extensive air showers are discussed. The 
main focus is put on the longitudinal development, in particular the energy 
flow, and on the lateral distribution of different air shower components. 
The intention of the paper is to provide a basic introduction to the subject 
rather than a comprehensive review. 
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1. Introduction 

Extensive air showers (EAS) are known for about 70 years to be cascades 
initiated by primary cosmic rays in the atmosphere. As such, they serve 
as connection to the highest particle energies nature is offering, especially 
at energies exceeding 10^^ eV where direct measurements of cosmic rays 
are hampered by the low primary flux. EAS can be viewed as tools for 
astroparticle physics. The main task is to reconstruct from the shower 
observables the parameters of the initial cosmic ray, which is in general 
straightforward for the particle direction. It is more a challenge for the 
energy and much more for the particle type of the primary. This is due 
mostly to shower fluctuations, an important EAS property given by the 
stochastic nature of the elementary interaction processes. Combined with 
the steep primary spectrum, the fluctuations impose a particular challenge 
for interpreting air shower data. Use can be made, however, of the fact that 
an EAS consists of different shower components. 
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Fig. 1. Schematic sketch of an EAS. 

The main body of this paper deals with illustrating, supported by EAS 
simulations, some characteristics of these shower components and the in- 
formation they contain about the primary particle. This might serve as 
introduction to the subject of extensive air showers and help to better un- 
derstand what air shower experiments are measuring, which information 
they try to infer from the data, and where limitations are given 01 • A more 
comprehensive introduction to EAS is given e.g. in jSj. 

A schematic sketch of an EAS is given in Figure ^ The incident particle 
(primary energy £"0 and primary type ^0) hits an air nucleus. Usually con- 
sidered as extreme illustrations of parimary particle types are protons and 
iron nuclei (although an analysis of EAS data is not restricted to this range) . 
The notation Aq anticipates the fact that EAS features of primary hadrons, 
which are known as dominant cosmic rays from direct measurements at 
energies below 10^^ eV, show some dependence on the nucleon number A 
of a nucleus, as explained later. In the first interaction, a number of sec- 
ondary hadronic particles are generated. Mostly by decay of the neutral 
pions, electromagnetic sub-cascades are initiated, and by decay of (part of) 
the charged pions, muons and neutrinos are produced. While neutrinos are 
hardly detectable in classic air shower experiments, muons might reach the 
observation level even from large altitudes. It is important to note that a 
competition of the charged pions between decay and interaction exists which 
depends on the pion energy and the local air density. Some information on 
the first interactions will be imprinted in this way in the muon component. 
The remainder of the hadronic secondaries continues to interact with air in 
subsequent collisions and to feed the other shower components. 

In ground arrays, the surviving particles are measured with particle de- 
tectors, and with appropriate optical telescopes Cherenkov and fluorescence 
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Fig. 2. Energy flow in EAS as obtained by CORSIKA shower simulations for an 
individual 10^^ eV primary proton event. The energy fractions stored in hadrons, 
electromagnetic particles, muons, and neutrinos are shown. The difference between 
their sum to the initial energy indicates the total amount of energy already released 
into air (shaded area). Upper graph in linear, lower graph in logarithmic scale. 



emission during the shower propagation can be observed. 



To characterize more quantitatively the interaction process, it is helpful 
to introduce the quantity inelasticity k as the energy fraction available for 
the production of secondary particles (or in other words, the initial collision 
energy reduced by the energy of the most energetic particle). As a rule 
of thumb, about is "lost" to the electromagnetic channel per hadronic 
interaction. For a mean value k ~ 0.6 typical for high-energy interactions 
this corresponds to ~ 20% per interaction that on average will give rise to 
subsequent electromagnetic cascading. Equipped with this knowledge, we 
can turn to the longitudinal shower development. 
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2. Longitudinal shower development 

2.1. Energy flow 

At first the general energy flow in an air shower is discussed in some de- 
tail. In Figure 121 the energy carried by different shower components during 
the cascade development are displayed for the example of a proton-induced 
event of primary energy 10^^ eV. These are results of detailed Monte Carlo 
air shower simulations obtained with the program package CORSIKA jH). 
For more information on simulation in particular at the highest energies, 
see for instance ^ ^'^d references given therein. 

Initially, all energy is concentrated in the primary proton. For this 
event, the first interaction occurs after traversing an atmospheric depth of 
~40 g cm~^. A significant energy fraction is transferred to the electromag- 
netic component, which is further fed in the subsequent shower process. 
One can notice an exponential decrease of the energy left in the hadrons. 
This can be understood in a simple picture of a constant inelasticity k and 
a constant mean free path for hadronic interactions. Based on the aver- 
age fraction of that is put into the electromagnetic channel per hadronic 
interaction, the energy remaining in the hadronic component at depth X is 

Eh{X) = Eo-{l-h)''/^'^ . (2.1) 

Adopting typical values used for modeling nucleon-air interactions at high 
energies of A/i ~ 55 g cm~^ and k ~ 0.6, the hadronic scale depth A/j of the 
exponential fall-off thus amounts to 

A,, = - — ^ , ~ 250 g cm-2 (2.2) 
|ln(l-iA;)| ^ ^ ^ 

in reasonable agreement with the results plotted for the detailed simulation. 

In this approximation, energy losses into the muon and neutrino channels 
were assumed to be small. In fact, as can be seen only at later stages 
of the shower development, an energy fraction of around 5% of the initial 
energy has accumulated in these components. The different character of the 
longitudinal curves of these two components with respect to hadrons and 
electromagnetic particles is evident. It is an important feature of muons to 
transport most of their (though small) energy fraction down to observation 
level. 

The electromagnetic particle channel is rising fast. The process to ex- 
tract energy from the electromagnetic particles is the energy loss due to 
ionization of air induced by the charged shower particles 0. Therefore, 
in Figure |31 the number of shower electrons (a term which usually includes 
positrons) is overlayed to the previous energy fiow plot. Since the specific 
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Fig. 3. Energy flow in EAS (left scale, compare Fig.|21) and cascade profiles N{X) or 
1 / a ■ dErei{X) / dX (right scale). The depth where the shower reaches its maximum 
number of particles is indicated by Xmax- 



ionization loss a of relativistic electrons is about constant (a ~ 2 MeV/(g 
cm~^)), the same curve in appropriate units also resembles the differential 
energy release dErei{X)/dX of the EAS as a whole, 

^{X)c.aN{X) . (2.3) 

Integrating dErei{X) / dX results in the energy fraction indicated by the 
shaded region; for instance, at shower maximum (labeled Xmax)^ already 
~ 50% of the initial energy has been released into the atmosphere. 

The maximum of energy stored in electromagnetic particles Xe/m — 
410 g cm~^ is reached well before the so-called shower maximum X^ax — 
730 g cm~^, i.e. the depth where the shower contains the largest electron 
multiplicity. This is due to the fact that at early cascade stages, a large 
energy fraction is carried by high-energy particles. Only gradually, the 
energy is transformed to newly created particles. The maximum of energy 
stored in electromagnetic particles Xg^m is expected at the development 
stage where the gain from the hadron channel equals the loss by energy 
release, 

-^(^e.™)-^(^e.™) . (2.4) 

We can roughly cross-check X^^im in our simplified approach. Using 
Eq. (1231) and (by differentiating Eq. jSH) -dEh{X)/dX ~ Eh{X)/Kh we 
obtain from Eq. (|2.4|) for X^im the condition 

Eh{Xelm)^aKhN{Xel^) . (2.5) 
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Adopting (rounded) values yields an expectation of -E'/i(410 g cm~^) ~ 2 • 
250 • 2 ■ 10^ MeV = 10% • Eq which is in reasonable agreement with the data. 

Equations H2.4|) and ()2.5() indicate the connection between energy flow 
(of hadronic to electromagnetic channel) and particle multiplication (in the 
electromagnetic channel). Most important for the observation are the pro- 
files dErei{X) / dX or N(X), respectively. These profiles are observable as 
fluorescence light and indicate how many particles reach the observation 
level. We therefore next investigate some characteristics of cascade curves 
formed by particle multiplication. 

2.2. Profile features 

Some of the main features of shower profiles can be nicely motivated 
within a simple toy model of particle cascades [H] , which complements to the 
previous toy model of energy flow. Let us suppose a particle with energy Eq 
that splits its energy equally into two particles after traveling a path-length 
A, and let this process be repeated by the secondaries. We then obtain a 
particle cascade which at a path-length X has evolved into N{X) = 2^1^ 
particles of equal energy E{X) = Eq/N{X). Let us further assume that 
particle multiplication stops when a certain energy limit E = Ei is reached. 
Then, the maximum number of particles N^ax is reached at this point X^ax^ 
and it is given by Nmax = N{Xmax) = Eq/Ei. The position of Xmax follows 

asX^a^ = A/ln2•ln(^o/^^)• 
Within this toy model, let us construct the more general case of an initial 
set {Aq,Eq/Aq) of ^0 particles, each with energy Eq/Aq. (The previous 
case is then realized for Aq = 1.) The particle number at maximum and the 
position of shower maximum are then given by 

Nmax{AQ, Eq/Aq) = NmaxiEo) OC Eq , (2.6) 

XmaxiAo, Eo/Ao) « ln(SoMo^O < X„,ax{Eo) . (2.7) 

If we identify the initial particle set (Aq, £"0/^0) as primary nucleus of mass 
number Aq (this is known as "superposition model"), the toy model predic- 
tions can thus be summarized as 

• Xmax increases proportional to the primary energy 

• Xmax is the same for all nuclei 

• Xmax increases with the logarithm of the primary energy 

• Xmax is smaller for the heavier nuclei (logarithmic dependence on Aq) 

• Xmax is the same for same Eq/Aq but different Eq 
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Fig. 4. Average depth of shower maximum as a function of primary energy [Jj. 



Despite the obvious simplicity (and its limitations) of this toy model, the 
findings of detailed shower simulations are quite well reproduced, as illus- 
trated in Figures \^ and [HI The average ^max increases with In Eq and is 
smaller for primary iron {Aq = 56) compared to primary proton, with a 
difference of ~ 80—100 g cm~'^ (Fig. [IJ. The primary iron Xmax for Eq fits 
well to the proton one of energy So/56 (same Eq/Aq). Nmax is quite sim- 
ilar for the different primaries (Fig. [S]) and about proportional to Eq (not 
shown) . 

Two other important EAS properties should be pointed out from Fig. |SJ 
Firstly, for fixed Eq different Xmax values translate into different parti- 
cle numbers on ground. For instance, the primary proton showers result 
on average in a larger number of ground particles compared to iron show- 
ers if the observation level is beyond the maximum. Thus, observation 
of Xmax or related quantities provide information on the primary particle 
type. Secondly, however, shower-to-shower fluctuations are visible, that lead 
to partly overlapping distributions of shower observables. As an example, 
the RMS (Xmax) of primary protons is ~ 60 g cm~^ and thus nearly as large 
as the difference between the average values of proton and iron. This limits 
an event-by-event assignement of a primary particle type, and composition 
analyses are usually performed with large event samples. However, it can 
also be seen that iron shower fluctuations are smaller (e.g. Y{MS>{Xmax) 
~ 20 g cm~^) than the proton ones. This might be used in turn to conclude 
about the primary composition by studying the very fluctuations of shower 
observables in a given event sample. 
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Fig. 5. Individual longitudinal shower profiles (vertical incidence). 

2.3. Shower electrons and muons 

The total number of shower particles regarded so far was dominated by 
the shower electrons. The shower muons, however, provide complementary 
information on the primary particle. As visible from the longitudinal dis- 
tributions in Figure (left), the muon particle number is decreasing only 
slowly after the maximum, in contrast to the shower electrons. Moreover, 
the total muon number also depends on the primary particle type. In the 
superposition model, the larger total muon content in iron showers might 
be qualitatively understandable: Due to the smaller energy per nucleon 
(Eq/Aq), the secondary pions are less energetic. This favours a pion decay 
as well as the fact that iron events develop at larger altitudes, where the air 
density is smaller. 

Shown in Figure IHl (right) are muon versus electron number for proton 
and iron induced events for different fixed primary energies. The size of each 
"potato" corresponds to the shower fiuctuation, while the separation indi- 
cates that a correlated measurement of ground particle numbers of shower 
electrons and muons allows conclusions on the primary mass. 

3. Lateral distribution 

Only the longitudinal development was discussed so far. Air showers 
have a lateral spread that also differs for the different shower components 
as well as for the various primary particles. It can be seen in Figure [7| (left) 
that the lateral distribution of shower muons on ground is flatter than the 
distribution of shower electrons. This is mostly due to the muon origin from 
larger altitudes compared to the more local production and fading of the 
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Fig. 6. (Left) Longitudinal profiles of typical primary proton and iron events for 
shower muons and electrons. {Right) Distribution of total muon and electron num- 
ber on ground for proton and iron induced showers of different primary energies [H| ■ 

electron component. In spite of the electrons being much more numerous 
than muons (about two orders of magnitude around shower maximum, see 
Fig. 121), at larger distance from the shower core the particle densities become 
comparable. 

Given the differences in the longitudinal development between primary 
proton and iron events, differences also in the lateral distributions might 
be expected. In Figure [3 (right), the ratio of particle densities of proton- 
induced to iron- induced events is displayed for shower muons and electrons. 
Both ratios decrease with increasing distance from the shower centre, in- 
dicating the flatter lateral distributions in case of the (on average) more 
developed primary iron showers. The larger muon content of iron-induced 
events is also visible (ratio < 1). The larger electron number on ground 
for primary proton showers can now be specified as larger ground particle 
density closer to the core, while at larger distances (which contribute less 
to the integrated, total particle number), the electron density in iron events 
is larger due to the flat lateral distribution. Thus, measuring local particle 
densities of the different shower components as a function of core distance 
provides additional information on the primary particle, which is utilized in 
large arrays of ground particle detectors. 

4. Discussion and conclusion 

Extensive air showers consist of different particle components which have 
different shower characteristics. This gives a handle to determine primary 
energy and (to some extent) the primary particle type. Extensions of this 
approach that were not discussed comprise the exploitation of different time 
structures of the shower front or comparing inclined events (where mostly 
muons survive the increased atmospheric path-length) to near-vertical ones. 
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Fig. 7. (Left) Lateral distribution on ground of shower muons and electrons. (Right) 
Ratio of lateral distributions on ground in proton to iron induced events, both for 
shower muons and electrons. 



Also in the hadronic shower component, information on the primary particle 
type is imprinted; but maybe even more important is the possibility to test 
high-energy hadronic interaction models used in the simulations. Shower 
simulations have grown to an important tool for EAS data reconstruction. 
Shower fluctuations motivate the development of Monte Carlo simulation 
techniques. In turn, these simulations might also be applied to find ob- 
servables that depend less strongly on the specific primary or on shower 
fluctuations, such as the signal at distances of 600—1000 m from the core 
for primary energy estimations with ground arrays. At the highest energies, 
further interaction features have to be considered in shower simulations, 
e.g. for primary photons the LPM effect and preshower formation in the 
geomagnetic field. 

In summary, even after many years of extensive exploration, air showers 
continue to be fascinating tools for astroparticle physics. 
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